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^ ! ABSTRACT 

D ■ 

C/} ■ The Microvariability and Oscillations of stars (MOST) satellite has detected 

\Q • low-amplitude light variations (Am ~1 mmag) in the Be star (3 CMi (B8Ve). The 

observations lasted 41 days and the variations have typical periods ~ 0.3 days. 
We demonstrate that the dominant frequencies are consistent with prograde high- 
order g-modes of m = — 1 excited by the Fe-bump of opacity in an intermediate- 
mass (~ 3.5M Q ) star with a nearly critical rotation period of 0.38 days. This 
Q\ ■ is the first detection of nonradial g-mode pulsations in a Be star later than B6 



leading to the possibility that pulsations are excited in all classical Be stars. 
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1. INTRODUCTION 

Be stars are non-supergiant, rapid rotators, which at times show emission in the Balmer 
and certain metallic lines. They are often called "classical Be stars" to distinguish them 
from other emission line stars such as Herbig Ae/Be stars and Algol systems (see Porter & 
Rivinius 2003, for a recent review on Be stars). The line emission is assumed to arise in a 
geometrically thin disk consisting of matter ejected from the central star. Such circumstellar 
disks have actually been imaged by long-baseline interferometry for several Be stars including 
P CMi (e.g. Tycner et al. 2005). Our target star, p CMi (HD 58715, HR 2845), is a bright 
but 'quiet' late-type Be star (B8Ve; V = 2.886, B — V = —0.072); no photometric variations 
were found by Pavlovski et al. (1997). The star shows only gradual and small long-term 
variations in Ha emission (e.g., Hanuschik et al. 1996) arising from a relatively small disk of 
~ 3 stellar radii (e.g. Tycner et al. 2005). 

Other Be stars, particularly of early-type, often show short-term photometric variations 
of less than a few days. The causes of these variations are assumed to be pulsation and 
rotational modulation (see Porter & Rivinius 2003, for a review). Many early-type Be stars 
exhibit line-profile variations indicating the presence of non-radial pulsations (e.g. Rivinius 
et al. 2003). Multi-periodic photometric variations are also the signature of non-radial pul- 
sations such as those found recently by MOST for ( Oph (09.5Ve; Walker et al. 2005a) 
and HD 163868 (B5Ve; Walker et al. 2005b). Until now, there has been no clear indication 
of pulsations for late-type Be stars. Baade (1989) searched for, but failed to detect any, 
line-profile variations in B8-B9.5 stars. There were very small photometric variations in a 
few late-type Be stars, but these were judged to be mono-periodic (e.g., \i Pic B9Ve, Balona 
et al 1992; ( Crv B8Vne, Barrera et al. 1991) and attributed to rotational modulation. 

It is of critical importance to know if nonradial pulsations are excited even in late-type 
Be stars because it would indicate that all Be stars may have nonradial pulsations which 
could play a critical role in mass ejection. Excitation of g-modes in late-type Be stars is 
likely because the cool boundary of the distribution of Be stars on the HR diagram (Zorec 
et al. 2005) roughly coincides with the cool boundary of the Slowly Pulsating B (SPB) stars 
(Pamyatnykh 1999), and because prograde g-modes can be excited even in a rapidly rotating 
star (Walker et al. 2005b). With this motivation we observed P CMi with MOST and indeed 
discovered multi-periodic light variations, albeit of low amplitude. 
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2. THE MOST PHOTOMETRY 

The MOST satellite was launched on June 2003 and the original mission is described 
by Walker et al. (2003). A 15/17.3 cm Rumak-Maksutov telescope feeds two CCDs, one for 
tracking and the other for science, through a single, custom, broadband filter (350 - 700 nm). 
Starlight from Primary Science Targets (V < 6) is projected onto the science CCD as a fixed 
(Fabry) image of the telescope pupil covering about 1500 pixels for high photometric stability 
and insensitivity to detector flatfield irregularities and the effect of particle irradiation on 
individual pixels. The experiment was designed to detect photometric variations with periods 
of minutes at micro- magnitude precision and does not rely on comparison stars or flat-fielding 
for the Fabry photometry. There is no direct connection to any photometric system. Tracking 
jitter was dramatically reduced by early 2004 to ~1 arcsec which led to significantly higher 
precision in the Fabry photometry. 

The observations received from the satellite were reduced by RK. Outlying data points 
generated by poor tracking or cosmic ray hits were removed. MOST suffers from parasitic 
light, mostly Earthshine, at certain orbital phases, with the amount and phase depending 
on the stellar coordinates, spacecraft roll and season of the year. Data are also recorded for 
Fabry images from seven of the eight lenses adjacent to the target Fabry lens in order to 
track the stray light background. These background signals were combined in a mean and 
subtracted from the target photometry. This also corrected for bias, dark and background 
signals and their variations. The reductions basically followed the scheme outlined earlier 
by Rucinski et al. (2004). 

MOST observed f3 CMi from December 16 2005 till January 26 2006 for a total of 41 
days. Observing time was split between (3 CMi and another field such that only 45% of each 
101 min orbit was available for (5 CMi. In addition to the 45% duty cycle over the 41 days, 
there was a single ~12 hour gap. Individual exposures were 8 sec taken every 20 sec (the 
sampling time). For the final frequency analysis mean magnitudes were calculated from data 
accumulated every 5 min. The RMS scatter of the individual 8-sec exposures is 0.300 mmag 
and 0.070 mmag for the the 5-min means. 

Figure 1 displays the full 41-day light curve for the 5-min means with time in heliocentric 
JD. A 3-day expanded portion of the light curve is shown in the second panel. The solid line 
in the lower panel is the reconstructed lightcurve obtained with the most significant peaks in 
the frequency analysis discussed below. The complete light curve can be downloaded from 
the MOST Public Data Archive K 



lr The MOST Public Data Archive is available at www.astro.ubc.ca/MOST. 
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Figure 2 displays the discrete Fourier transform (DFT) of the MOST light curve for 
/3 CMi. The frequency analysis was performed by CC using his CAPER software (Walker et 
al. 2005b; Saio et al. 2006). Frequencies with amplitudes > 0.04 mmag (3.5 times the noise 
level averaged over the spectrum from to 15 cycles day -1 ) are identified in Figure 2. The 
plot includes the spectral window function and shows the residual amplitude spectrum after 
removal of the 20 highest peaks. The peak at 1 cycles day -1 is an artifact of the Earth's 
rotation and the satellite's Sun-synchronous orbit, and their modulation of parasitic light. 
The peak at 15.2 cycles day -1 is an artifact related to MOSTs orbital frequency. 

There are two dominant features in the DFT, with frequencies of 3.257 ± 0.001 cy- 
cles day -1 (S/N ~ 8.5) and 3.282 ± 0.004 cycles day -1 (S/N ~ 3.5). The spacing between 
these two frequencies is 0.025 ± 0.004 cycles day -1 . The high S/N values for these fre- 
quencies and detailed comparison with the background recorded in the MOST Fabry field 
demonstrates that these two signals are intrinsic to f3 CMi. 

There are also nearby frequencies at 3.135±0.006 cycles day -1 (S/N ~ 2.9) and 3.380± 
0.002 cycles day -1 (S/N ~ 2.8). However, given their relatively low signal-to-noise, we do 
not consider these to be definitive detections in this data set. There is also evidence for 
significant stellar variability at frequencies below 3 cycles day -1 but its origin is unclear. 

We have performed the frequency analysis on the data, both with the long-term trends 
removed and left in the time series. The frequency identifications reported above are unaf- 
fected, which is consistent with the very clean spectral window of the MOST photometry and 
the expected lack of power leakage from low to higher frequencies in the Fourier spectrum 
of these data. 

Error bars for the fit parameters are determined by a bootstrap process (e.g., Wall & 
Jenkins 2003). A large number of trial light curves (10000 in this case) are generated by 
randomly sampling the (3 CMi light curve. The fitting procedure is repeated for each new 
light curve resulting in a distribution for each fit parameter. The error bars are 1 standard 
deviation in each of the parameter distributions (see for example Saio et al 2006). Note 
that these frequencies are well determined although the amplitudes of the main peaks have 
moderate uncertainties. 



3. MODELS 

To understand the cause of the periodic variations detected by MOST we have per- 
formed a pulsational stability analysis for rapidly rotating main-sequence star models of 3.5 
and 3.6M , whose evolutionary tracks pass close to the position of (3 CMi (Fig. 3). The 
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evolutionary models were computed in the same way as in Walker et al. (2005b), and the 
stability analysis was performed based on the method of Lee & Baraffe (1995). An initial 
chemical composition of (X, Z) = (0.70, 0.02) is adopted. The opacity was obtained from 
OPAL tables (Iglesias & Rogers 1996). The centrifugal force is approximately included in 
the spherical symmetric hydrostatic equilibrium as 

*=-»p+?«". (i) 

where P is the pressure, r is the distance from the center, g is the local gravitational accel- 
eration, p is the matter density, and Q is the angular frequency of rotation. The rotation 
frequency is assumed to be constant throughout the stellar interior and evolution. Evolu- 
tionary tracks calculated with rotation frequencies of 0.02 mHz (1.73 cycles day -1 ) and 0.03 
mHz (2.59 cycles day -1 ) are shown in Figure 3. These tracks stop where the centrifugal force 
exceeds the gravitational force in equation (1) at the stellar surface. A track not including 
rotational effects has also been calculated and is also shown in Figure 3. 

Figure 3 also shows an approximate position of (3 CM i with error bars, which is estimated 
as follows. Zorec et al. (2005) obtained (logT eff , logg) = (4.070,3.88) for j3 CMi as apparent 
parameters and converted them to mean values (logT c ff, \og~g) = (4.081,3.94) averaged over 
the whole stellar surface. These parameters yield logL/L = 2.33. On the other hand, 
Fabregat & Reglero (1990) obtained an absolute magnitude of M v = 0.01 from the data 
of uvbyP photometry corrected for circumstellar emission. Combining this value with a 
bolometric correction for B8V of —0.7 (Flower 1977) yields logL/L Q = 2.18. Furthermore, 
combining Fabregat & Reglero's V magnitude of 3.02 (corrected for circumstellar emission) 
with the Hipparcos parallax (Perryman et al. 1997) and the above bolometric correction, we 
obtain log L/L Q = 2.41. From these estimates we have adopted log L/L Q = 2.29 ± 0.11 for 
the luminosity of f3 CMi. For the effective temperature, Flower (1977) gives logT e ff = 4.086 
for the spectral type B8V. Combining this value and Zorec et al.'s value, we have adopted 
logT cfr = 4.081 ±0.005. 

(3 CMi is located close to the cool edge of the SPB instability region of the HR dia- 
gram (Pamyatnykh 1999). Our models confirm this observation (see Figure 4). Along our 
evolutionary tracks without rotation the cooler boundary of the instability region appears 
at logT efr « 4.03 for I = 1, ta 4.06 for I = 2, and « 4.09 for I = 3. The excited modes 
in our 3.5M , non-rotating models, are high-order g-modes excited by the Fe-opacity bump 
(Gautschy & Saio 1993; Dziembowski et al. 1993). 

Since (3 CMi is a rapid rotator, we have to include the rotation effects in the stability 
analysis of g-modes. Nonradial pulsations of a rapidly-rotating star cannot be represented 
by a single spherical harmonic and coupled with toroidal velocity fields. In our analysis the 
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angular dependence of pulsation amplitude is expanded using eight spherical harmonics Y™ 
for a given azimuthal order m (Y™ for toroidal velocity field), where lj = \m\+2j (lj = lj + 1) 
for even modes and lj = \m\ + 2j + 1 (lj = lj — 1) for odd modes with j = 0, 1, . . . , 7 (see 
e.g., Unno et al 1989, for general discussions on nonradial pulsations of rotating stars). We 
adopt the convention that a negative m represents a prograde mode (in the co-rotating frame) 
with respect to the stellar rotation. Even (odd) modes are symmetric (anti-symmetric) with 
respect to the equatorial plane. We designate the angular-dependence type of a mode by a 
set of (m,£), in which £ is defined as the value of the largest-amplitude component. 

Figure 5 shows frequencies (in the observers' frame) of low surface order (\m\ < 2 and 
£ < 2) modes excited in evolutionary models of 3.5M Q rotating at frequencies of 0.02 mHz 
and 0.03 mHz as a function of the effective temperature. Comparing this figure with Figure 4, 
we find that the cooler boundary of instability is bluer than in the non-rotating case. Fewer 
modes are excited in more rapidly rotating models. This implies rapid rotation tends to 
stabilize g-mode pulsations. 

The pulsation frequency in the observers' frame is given as 

|i/(crot) - mi/ r ot|, 

where z/(crot) is the pulsation frequency in the frame co-rotating at frequency v xoXi . Therefore, 
pulsation frequencies in the observer's frame for high-order g-modes and r-modes (for which 
i/(crot) <C |m|z/ rot ) tend to group depending on the m value just as seen in Figure 5. All the 
excited prograde modes (m < 0) seen in Figure 5 are sectoral (£ = \m\) g-modes, while all 
the excited retrograde (m > 0) modes are r-modes (in which toroidal motion is dominant) 
in those rapidly rotating models. No retrograde g-modes are excited in this figure. We note 
that, as shown in Appendix A, the stability of pulsation modes differs considerably if we use 
the so-called traditional approximation in which the horizontal component of the angular 
velocity of rotation is neglected. 

Under the assumption of u rot m 0.02 mHz (1.73 cycles day" 1 ), the periodic light varia- 
tions of (3 CMi are approximately reproduced by r-modes of (m,£) = (2,2). On the other 
hand, the observed variations are identified as sectoral prograde g-modes of (m,£) = (—1, 1) 
if we assume z/ rot ps 0.03 mHz (2.59 cycles day" 1 ). In the former case another three groups 
of frequencies are expected to be observed in the estimated logT eff range of (3 CMi, while 
in the latter case only one group of frequencies should be observed (Fig.5). Since MOST 
has detected only one group of frequencies around 3.26 cycles day -1 in (3 CMi (§2), we 
conclude that a model rotating at u rot 0.03 mHz (2.59 cycles day" 1 ) is more appropri- 
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ate. 2 For this rotation rate pulsations excited at logT cff = 4.083 are prograde sectoral 
high-radial order (19th-33rd) g- modes. The frequencies in the co-rotating frame range over 
0.586 cycles day -1 > z/(crot) > 0.334 cycles day -1 . Frequency spacings between two adja- 
cent modes lie in a range from 0.035 cycles day -1 to 0.009 cycles day -1 . Since the frequency 
resolution expected from the data is approximately 1/(41 days) = 0.024 cycles day -1 , only 
a few frequencies are expected to be resolved just as the DFT indicates (Fig. 2). 

The top panels of Figure 5 show the rotation angular frequency Q normalized by the 
Keplerian frequency (GM/R 3 ) 1 ^ 2 . For a given z/ rot , f2 increases as the model evolves because 
the stellar radius R increases. Since the equatorial radius is larger than the mean radius, 
the rotation velocity becomes critical at the equator with Q < 1; probably at Q ~ 0.8. 
Figure 5 indicates that 3.5M models with z/ rot = 0.03 mHz rotate nearly critically at the 
effective temperature of (3 CMi. (For 3.6M models with z/ rot = 0.03 mHz, f2 > 0.8, not 
appropriate for (3 CMi.) This supports Townsend et al's (2004) claim that classical Be stars 
may be rotating much closer to their critical velocities than is generally supposed. It is also 
consistent with the recent statistical property found by Cranmer (2005) that the coolest Be 
stars rotate nearly critically. 

4. DISCUSSION 

We have found nonradial g-mode pulsations excited in the late-type Be star (3 CMi 
(B8Ve). This confirms that f3 CMi belongs to the new class of SPBe stars, like HD163868 
(Walker et al. 2005b). We have identified pulsation modes as prograde-sectoral g- modes of 
m — — 1, based on the comparison between the observed frequencies and theoretical ones 
excited by the Fe bump of opacity in rapidly rotating models. The fact that only one group of 
frequencies is detected by MOST is consistent with our prediction of the pulsational stability 
if (3 CMi rotates nearly critically. This supports the recent claims (Townsend, et al. 2004; 
Cranmer 2005) that classical Be stars rotate more rapidly than previously thought (70-80% 
of the critical rate; e.g. Porter 1996; Chauville et al. 2001). The observed amplitudes are 
small (< 1 mmag). Since (3 CMi lies close to the red edge of the instability region on the 
HR diagram, the pulsational instability is not strong and hence the pulsation amplitudes are 
expected to be small. 



2 The slower rotation model might appear to be saved if the low-amplitude frequencies around 2 cy- 
cles day -1 in Fig. 2 were assigned to g-modes of (m,£) = (—1,1) and the main frequencies around 3.26 
cycles day -1 to r-modes of (m,£) = (2,2). We consider this highly unlikely, however, since r-modes are 
not expected to yield luminosity variations much larger than the sectoral g-modes of (m,£) — (—1, 1), and 
because no peaks around 4 cycles day -1 corresponding to g-modes of (m,£) = (—2,2) are detected. 
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Our discovery of nonradial pulsations in (3 CMi suggests the possibility that nonradial 
pulsations are involved in all rapidly rotating Be stars. In some early-type Be stars mass 
ejections are observed to occur at certain pulsation phases (e.g. \x Cen: Porter & Rivinius 
2003) suggesting that pulsations trigger mass ejections. It is tempting to speculate that 
nonradial pulsations may play a crucial role in mass ejection in all Be stars. Osaki (1986) 
argued that prograde nonradial pulsations can transport angular momentum toward the 
surface and spin up the equatorial region critically leading to a mass ejection. Owocki 
(2005) showed numerically that introducing velocity variations corresponding to a prograde 
(m = —2) nonradial pulsation in a nearly critically rotating stellar atmosphere actually 
yields a circumstellar Keplerian disk. 

Our findings in this paper and Walker et al. (2005b), that prograde sectoral g- modes are 
selectively excited in rapidly rotating Be stars, support the above scenario for mass ejection. 
Line-profile-variation (LPV) analyses for early type Be stars, however, tend to indicate the 
presence of retrograde modes (Rivinius et al. 2003) rather than prograde modes, contrary to 
our model findings. 

Rivinius et al. (1998) obtained four periods around 0.503 days and two period around 
0.28 days from LPVs in \x Cen, one of the best studied and most active Be stars. [These 
LPVs of \i Cen were modeled as retrograde (m = 2, and 3) modes by Rivinius et al. (2001).] 
Frequency spacings among each group are about 0.034 cycles day -1 or less, which are com- 
parable with the frequency difference 0.025 cycles day -1 between the two main frequencies of 
/3 CMi (§2) and with typical frequency spacings ~ 0.03 cycles day -1 of HD 163868 (Walker 
et al. 2005b). These frequency spacings are comparable with those of g-modes, although 
they are close to the limit of the frequency resolution. Longer observations by MOST or 
COROT would make g-mode seismology of SPBe stars possible. Also, simultaneous pho- 
tometric and spectroscopic observations of ji Cen and other Be stars showing conspicuous 
LPVs and mass-ejection events would be important to understand the connection between 
nonradial pulsations and mass ejections in Be stars. 
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research has made use of the SIMBAD database, operated at CDS, Strasbourg, France. The 
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(Quebec). R.K. is supported by the Canadian Space Agency. W.W.W. is supported by the 
Austrian Space Agency and the Austrian Science Fund (P17580-N2). H.S. is supported by 
the 21st Century COE programme of MEXT, Japan. 
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A. Mode stability under the traditional approximation 

The pulsation analysis for low-frequency modes of a rapidly rotating star is greatly sim- 
plified by the use of the so-called traditional approximation, in which the horizontal com- 
ponent of the angular velocity of rotation (fl sin 9 with 9 being the co-latitude) is neglected. 
Under the traditional approximation, eigenfrequencies can be obtained using the same equa- 
tions as those for nonradial pulsations of non-rotating stars except replacing £(£ + 1) with 
the parameter (see e.g., Townsend 2005; Lee 2006). The value of deviates from 
£(£ + 1) as the ratio of the rotation to the pulsation frequency (in the co-rotating frame) 
increases depending on to and £ (e.g., Lee & Saio 1997). 

The traditional approximation is justified dynamically for low-frequency pulsations, but 
it is not clear whether the approximation is reasonable for the nonadiabatic stability analysis. 
To see the effect, we have performed a stability analysis using the traditional approximation 
for 3.5-M models rotating at a frequency of 0.03 mHz. The results are shown in Figure 6. 
Comparing this figure with Figure 5, we see that under the traditional approximation, larger 
numbers of prograde sectoral g-modes are excited, and that tesseral g-modes of (to, £) = (0, 1) 
and (—1,2) are also excited. These tesseral modes are all damped without the traditional 
approximation in these models. The cause of the difference is attributed to damping due 
to mode-couplings (Lee 2001) which are absent under the traditional approximation. The 
damping effect is strongest for retrograde g-modes and tesseral (£ > \m\) g-modes. 

We note that in the models shown in Figure 6 no retrograde g-modes are excited even 
under the traditional approximation. This is understood as follows. Retrograde g-modes have 
large values of (Lee & Saio 1997), which correspond to modes of large latitudinal degree £ 
in a non-rotating model. Figure 4 for nonrotating models indicates that the cooler boundary 
of instability region is hotter for larger £. Therefore, we understand that the temperature 
range for instability of retrograde g-modes is hotter than that shown in Figure 6. 
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Fig. 1. — a) MOST 41-day light curve of (3 CMi. Each point is a 5-min mean and the 
point-to-point RMS scatter is 0.070 mmag. b) An expanded three day portion of the light 
curve in (a). The solid line is the fit of the 20 highest peaks from the frequency analysis of 
the full light curve. 
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Fig. 2. — a) (top panel) shows the Fourier amplitude spectrum of the light curve in Figure 
1. The points with bars are frequencies returned by CAPER having the largest amplitudes. 
The window function is shown for the strongest peak. Residuals after removing the 20 largest 
peaks are shown at the bottom of the panel (the sign is inverted for better visibility). Small 
dots indicate the frequency and amplitude obtained at each bootstrap (e.g., Wall & Jenkins 
2003), in which the frequency analysis by CAPER was repeated many times for data sets 
consisting of randomly chosen data points (see Saio et al. 2006, for details). Dash-dotted lines 
indicate a level of 3.5 times the noise averaged over the spectrum from to 15 cycles day -1 , 
b) (bottom panel) shows a zoomed region around the most prominent peak. Note that all 
identified frequencies are resolved. 
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Fig. 3. — Evolutionary tracks of 3.5M and 3.6M with various rotation frequencies. Circles 
indicate positions of zero-age-main sequence models. Tic marks along each evolutionary track 
indicate positions separated by an interval of 5 x 10 7 yr. The square with error bars indicates 
the approximate position of j3 CMi. 
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Fig. 4. — Frequencies of g-modes excited in non-rotating models along the 3.5M main- 
sequence evolution. The evolutionary track is shown in Fig. 3. 
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Fig. 5. — Frequencies in the observer's frame of excited pulsation modes (bottom panels) and 
the angular velocity of rotation normalized by (GM/ R 3 ) 1 ^ 2 (top panels) as functions of effec- 
tive temperature. Horizontal dotted lines indicate observed frequencies 3.257 cycles day -1 
and 3.282 cycles day -1 for (3 CMi. Vertical dash-dotted lines indicate the estimated range for 
the effective temperature of (3 CMi. Rotational frequencies are assumed to be 0.02 mHz (left- 
hand panel) and 0.03 mHz (right-hand panel) through the main-sequence evolution. Symbols 
indicate type of modes: filled circles (•) for m = 0,£ = 1; triangles (A) for m = —1,1 = 1; 
squares (□) for m = —2,£ = 2; inverted triangles (▼) for m = l,£ — 2; diamond (o) for 
m = 2,£ = 2. Only r-modes are excited among retrograde modes (m > 0). Other excited 
modes are sectoral prograde g-modes. (A few axisymmetric g-modes are also excited in the 
hottest models). 
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Fig. 6. — The same as the right-bottom panel of Figure 5 but showing results obtained under 
the traditional approximation. Filled triangles (A) are for g-modes of (m,£) = (—1,2), and 
other symbols are the same as in Figure 5. Compared to Figure 5 more sectoral g-modes of 
m = —1, and —2 are excited. Tesseral g-modes of (m,£) = (0, 1) (around 1.7 cycles day -1 ) 
and of (m,£) = (—1,2) (around 4.7 cycles day -1 ) are excited only under the traditional 
approximation in these models. 



